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Dynamical friction is the process responsible for matter transport toward the inner region 
of galaxies in form of massive objects, like intermediate mass black holes, globular clusters 
and small satellite galaxies. While very bright galaxies show an almost flat luminosity 
profile in the inner region, fainter ones have, usually, a peaked, cuspidal, profile toward 
the center. This makes unreliable, in these cases, the use of the classic Chandrasekhar's 
formula for dynamical friction in its local approximation. Using both A'^-body simulations 
and a semi analytical approach, we have obtained reliable results for the orbital decay of 
massive objects in cuspidal galaxies. A relevant result is that of a shallower dependence 
of dynamical friction braking on the satellite mass than in the usual Chandrasekhar's 
local expression, at least in a range of large satellite masses. 
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1. The dynamical friction process 

The deceleration on a test mass M at position identified by the radius vector r 
exerted by the 'sea' of particles of mass m where it is moving in, is given by the 
scattering integral over the whole phase space occupied by the particles' : 
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where b is the impact parameter of the M~m two-body encounter. To evaluate this 
complicated integral along the M test particle path, the so called local approxima- 
tion formula for dynamical friction (df) is widely used'^ : 

= -A7:'G\m + M)p{r)F{r,VM)\nA^, (2) 

at Vjyj 

where F{r,VM) is the fraction of scatterers slower than the test particle and 
In A is the Coulomb's logarithm. This expression gives good results also in central 
galactic cores, but fails when the test particles moves towards the center of a galaxy 
host whose density profile shows a central cusp, as demonstrated by.'^ Actually, 
the central density divergence implies a significant overestimate of the df effect on 
objects transiting regions close to the center. To avoid this problem, Arca-Sedda 
and Capuzzo-Dolcetta,'^ in the hypothesis of spherical symmetry for the host galaxy, 
give an estimate of df by mean of the interpolation formula 
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where the cen and loc decelerations are obtained via the numerical evaluation of 
the integral in Eq. [T] and via the usual local approximation, respectively. In Eq|3] 
the weight function, p{r) is assumed as p{r) — exp{—r/rc), where Vc is a proper 
scale length. 

A more extended and deeper discussion of the treatment of dynamical friction 
in cuspy galaxies and of the validation of the method described above is found in.'^ 
In the following we just give a short summary of main aspects of the problem and 
some relevant results. 

2. Massive objects decay in cuspy elliptical galaxies 

To test the approach presented above, we performed a set of iV-body simulations 
of the motion of massive objects in galaxy models assumed as'^ spherical density 
profiles, characterized by a length (r^) and a mass (Mg) scale, and diverging as 
p(r) cx r"'*' in the inner region. 

The simulations were done at varying: i) the satellite mass, 
M/Ms (5 X 10"'', 10"^ 5 X 10^3), ii) the 7 value, 7 = (1/2, 1, 3/2), iii) the initial 
position (0.2 < ro/r^ < 2), and, iv) the type of orbit. The TV-body simulations were 
performed by using HiGPUs, a 6*'* order Hermite integrator running on parallel 
CPU+GPU systems, developed by.^ On the other side, the one body problem of 
motion of the satellite of mass M in the Dehnen's potential with the semianalytic 
treatment of dynamical friction described in Sect. [T] was integrated by mean of a 
Qth Qj.j(,j. Kungc-Kutta-Nystrom method. 

The comparison of the direct A^-body and semianalytic results over a set of test 
cases show that the difference in evaluating the decay time is less than 1% for radial 
orbits and below 10% for circular orbits. This indicates as possible and reliable to 
simulate the satellite orbital decay by mean of the semyanalitic approach avoiding 
the use of direct iV-body simulations, which would require a significant amount of 
computational time when covering large inetervals of the allowed initial conditions 
and parameters. 

The most important output is the decay time (tdf), that is the time the object of 
mass M requires to reduce its motion to a small oscillation around the center of the 
host galaxy. The local approximation formula leads to an inverse mass dependence 
for the decay time, tdf cx M~^. 

The analysis of our results, however, shows that the decay time has a somewhat 
different scaling with the satellite mass: 

tdf oc M-°-«^ (4) 

3. Results and applications 

Just to show some preliminary results, we show in FiglT] the df time for circular 
and radial orbits starting at the same position, ro/r^ = 0.2 — 2, for several masses 
decaying in a galaxy model of given mass Ms- The satellite mass range covers the 
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GC mass range up to dwarf spheroidals (to take into account minor mergers). The 
results shown in Fig. [1] are easily scaled to physical time units given the relation 
expressing crossing time r — rl^"^ [GMg)^^^'^ ■ 

As an example showing the relevance of dynamical friction in cupy galaxies, let 
us consider a somewhat standard galaxy as having Mg = 5 x 10^^ M0 and — 1 
Kpc. This leads to the age of the Universe (assumed 13 Gyr) to be, in logarithm, 
4.3 in the units of Fig. [l] i.e. in the uppermost part of that Figure. 




Fig. 1. Dynamical friction decay time (in units of crossing time t) as function of the satellite 
mass. The three different areas in the plot are delimited by a lower boundary (radial orbits) and 
an upper (circular orbits) and correspond to test masses start moving at the initial distance, ro 
in units of r^, as labelled at top right. All values refer to the Dehnen's 7=1 model. 



